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F-18  ESCAPE  SYSTEM  COMPUTER  SUPPORT 
FINAL  REPORT 

1.0  INTRODUCTION 

The  following  report  documents  the  work  completed  by  Computer  Sciences 
Corporation  (CSC)  under  Task  Order  Number  22,  Contract  N62269-75-C-0001  for 
the  Crew  Systems  Department  at  the  Naval  Air  Development  Center  (NADC), 

Warminster,  Pennsylvania.  This  report  documents  the  modifications  and  additions 
that  were  made  to  an  existing  emergency  egress,  pilot  ejection  from  an  aircraft, 
simulation  model.  The  original  model,  nicknamed  ICARUS  by  the  originators,  was 
developed  and  implemented  in  FORTRAN  on  a CDC  6600  computer  by  the  Aeroballlstics 
Division  of  the  Naval  Weapons  Laboratory  (NWL)  located  in  Dahlgren,  Virginia. 

Documentation  of  the  original  version  of  ICARUS  is  covered  under  NWL 
Technical  Report  Number  TR-3098,  dated  February,  1974  and  titled:  'Air  Crew 
Automated  Escape  System  Simulation  Model'  and  NWL  Technical  Note  TN-K- 7-7/74 
dated  February,  1974  and  titled:  'Program  Maintenance  Manual  for  ICARUS  A 
Computer  Program  to  Simulate  Escape  Systems. ' 

The  original  model  was  developed  to  simulate  the  ejection  seat  system 
of  the  AV-8A  (Harrier)  aircraft.  The  ejection  system  was  designed  and  manufactured 
by  the  STENGEL  Corporation.  The  model  was  constructed  and  programmed  in  a 
modular  form,  thus  enabling  the  current  or  future  users  of  the  model  to  modify  or 
insert  new  modules  into  the  model  in  order  to  adapt  the  model  to  hardware  changes 
to  the  original  ejection  system  or  different  ejection  systems. 

The  basic  structure  of  the  model  can  be  considered  as  follows: 

• Model  Control  Section  - This  section  controls  the  overall  model 
data  initialization,  the  individual  case  study  data  initialization, 
requests  for  aircraft  and  the  seat  ejection  normal  trajectory 
computations,  trajectory  data  output,  adjustments  to  the  trajectory 
time  step  (deltat)  and  case  or  problem  termination. 
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• Poeltlonal  or  Time  Event  AcMvation/Deactlvatton  Section  - This 

section,  colled  by  the  above  section  after  each  step  of  the  trajectory 

computation,  keeps  track  of  each  of  the  position  and  time  event's 

I status  during  the  trajectory.  When  an  event  occurrence  has  been 

i 

detected  by  this  section,  the  section  initiates  a request  to  have  a 
message  outputted  indicating  that  the  particular  event  has  taken 
place,  the  time  of  occurrence  to  be  computed  and  outputted  along 
with  the  appropriate  trajectory  data  computed  at  that  time.  This 
section  is  also  responsible  for  the  dynamic  scheduling  of  the  var- 
ious equations  of  motion  to  be  evaluated  at  each  phase  of  the  tra- 
I jectory  computation. 

i 

I The  first  request  from  the  NADC  Crew  Systems  Department  was  that 

\ CSC  change  the  output  format  of  the  model  (ICARUS)  from  a paragraph  to  a tab- 

^ ular  format.  A series  of  new  reports  were  designed  displaying  the  trajectory 

, data  in  tabular  form  (up  to  nine  components  per  page)  with  appropriate  page  and 

title  headings.  The  input  routine  was  modified  to  read-in,  immediately  after  the 
[ case  label  card,  an  output  display  selection  card.  The  input  routine  was  further 

’ modified  to  interpret  the  data  read-in  from  this  card  and  to  set  program  transfer 

switches  within  the  program  to  cause  the  printing  or  inhabiting  of  the  printing,  of 
the  data  or  reports  selected  by  the  user  via  this  card.  This  card  enabled  the  user 
• of  the  model  to  select  either  the  original,  new  or  both  report  formats  to  be  gen- 

erated and  displayed.  If  the  new  report  format  mode  of  output  data  display  is 
selected  by  the  user,  the  ability  to  select  specific  reports  contained  within  the 
mode  was  also  provided  via  this  card.  Simultaneously,  the  original  input  routine 
was  segmented  into  smaller  routines  and  the  user  of  the  model  was  given  the 
option  (also  via  the  above  card)  of  selecting  which  input  data  segment  was  to  be 
displayed  and  whether  the  recovery  chute  data  and  plot  data  is  to  be  outputted. 

CSC  was  then  asked  to  modify  ICARUS  to  compute  and  display  the 
pilot's  dynamic  response  index  (DRI),  a measure  of  the  pilot's  spinal  compression. 

A spring-damper  system  is  used  as  a model  of  the  pilot's  spine  during  the  ejection 
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phase.  The  DRI  is  computed  as  a function  of  the  vertical  acceleration  imposed 
on  the  man  during  this  phase.  The  computation  was  implemented  in  the  model 
and  a new  report  was  designed  to  display  this  data  at  each  time  step  of  the  ejec- 
tion phase.  This  accomplished,  CSC  was  then  asked  by  Crew  Systems  to  modify 
the  model  to  take  into  account  the  displacement  of  the  pilot's  center  of  gravity  dur- 
ing the  ejection  phase  (from  catapult  Ignition  to  seat/man  separation).  The  dis- 
placement of  the  pilot's  center  of  gravity  was  modelled  by  Crew  Systems  using 
spring  damper  equations  (one  for  each  direction  of  motion).  These  equations 
were  inserted  into  the  model,  evaluated  and  the  appropriate  vectors  were  up- 
dated to  reflect  the  displacement  of  the  pilot's  center  of  gravity  from  its  original 
(static)  position.  A new  output  report  displaying  the  displacement,  acceleration  and 
force  components  in  each  direction  of  motion,  was  designed  and  implemented  in 
the  model. 

The  preceeding  modifications  to  the  model  are  considered  common  to 
all  ejection  systems.  The  following  modifications  have  been  implemented  in  the 
model  in  order  to  adapt  and  use  the  model  to  simulate  a different  ejection  system 
namely,  the  Martin-Baker. 

The  Martin-Baker  catapult  was  modelled  by  Crew  Systems  using  spring- 
damper  and  tube  extension  equations.  The  rail  Interaction  with  the  seat  slippers 
was  modelled  using  spring  equations  to  simulate  and  determine  the  friction  and 
normal  forces  and  moments  generated  along  each  axis  as  the  slippers  move  up 
the  rails  during  the  catapult  phase.  The  forces  and  moments  acting  on  the  seat/ 
man  combination  due  to  the  bending  of  the  tube  while  extending  were  modelled 
using  tube  bending  equations.  The  above  equations  were  interfaced  and  implemented 
in  the  model.  Additional  reports  displaying  the  slipper  and  tube  bending  position, 
force  and  moment  data  were  designed  and  implemented  in  ICARUS. 

The  model  was  then  modified  to  simulate  the  Martin-Baker  seat  back 
rocket,  drogue  and  recovery  chute  sequencing  logic. 
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Three  new  decimal  data  input  sections  were  added  to  the  model  allow- 
ing the  user  of  the  model  to  input,  via  the  model's  normal  input  data  procedures, 
the  parameters  to  the  dynamic  center  of  gravity,  the  Martin-Baker  cat:^>ult  and 
phasing  equations. 


The  following  sections  (not  necessarily  in  chronological  order  of  im- 
plementation) describe  in  detail  the  mathematical  modelling,  timing  and  seq- 
uencing logic,  the  input/output  changes  and  additions  to  the  model,  a description 
of  the  new  routines  that  were  implemented  and  interfaced  with  the  rest  of  the 
model,  macro  and  detailed  flow  diagrams  of  the  model,  a definition  of  the  new 
input  variables  and  how  to  enter  them,  giiidelines  as  to  how  to  go  about  implement- 
ing the  new  input  variables,  and  finally  how  to  run  the  model. 


2.0  MATHEMATICAL  MODELLING  CHANGES  AND  ADDITIONS 

2.1  Dynamic  Response  Index  (DRI) 

A measure  of  the  compressibility  of  the  pilot's  spine  during  the  ejec- 
tion phase  of  the  trajectory  has  been  implemented  and  interfaced  in  ICARUS. 

The  pilot's  spine  is  modelled  as  a mechanical  system,  a spring-damjaer.  The 
index  of  compressibility  (DRI)  of  the  man's  spine  is  determined  at  each  time 
step  of  the  phase  as  a function  of  the  vertical  acceleration  (in  G's)  imposed  on 
the  seat/man  combination  and  the  resulting  displacement  value  obtained  from  the 
spring  equation  during  this  phase.  Subroutine  DERIV  has  been  modified  to  call 
a new  subroutine  titled  'DRICALC  during  the  ejection  phase.  Subroutine  DRICALC 
evaluates  the  spring-damper  equation  using  the  velocity  and  position  data  obtained 
after  each  numerical  integration  iteration.  The  dynamic  response  index  is  then 
computed  using  the  displacement  value  obtained  from  the  integration. 

Reference  Appendix  A for  a mathematical  presentation  of  the  equations 
used  to  determine  the  DRI. 
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2.2 


Dynamic  Center  of  Gravity 


The  orignlal  model  was  modified  to  take  Into  account  the  displacement 
of  the  pilot's  center  of  gravity  during  the  ejection  phase  of  the  trajectory,  from 
catapult  Ignition  to  seat/man  separation.  The  displacement  of  the  man's  center 
of  gravity  Is  simulated  during  this  phase  by  three  independent  second-order 
differential  equations  of  motion  (spring-damper  equations)  one  for  each  direction 
of  motion.  The  equations  are  numerically  integrated  at  each  time  step  as  a func- 
tion of  the  resultant  axial  acceleration  Imposed  on  the  seat/man  combination  - 
yielding  the  corresponding  velocity  and  displacement  components  of  the  man's 
center  of  gravity  from  its  static  position.  The  seat/man  center  of  gravity  la  then 
recomputed  and  the  appropriate  vectors  from/to  this  new  position  are  then  updated 
in  the  seat/man  coordinate  system. 

The  user  of  the  model  has  been  given  the  option  of  forcing  the  model  to 
compute  the  displacement  of  the  man's  center  of  gravity  or  ignoring  this  displace- 
ment as  was  originally  done  In  the  model.  This  option  is  controlled  by  the  user 
of  the  model  setting  the  fifth  word  of  the  integer  input  data  array  to  an  integer  one 
if  the  displacement  of  the  man's  center  of  gravity  from  its  static  position  is  to  be 
computed,  a zero  or  blank  entry  will  force  the  model  to  by-pass  this  computation. 

Subroutine  DERIV  has  been  modified  to  call  a new  subroutine  titled 
'DYNCG'  during  the  ejection  phase.  Subroutine  DYNCG  calls  another  new  sub- 
routine titled  'KUTTA'  which  performs  the  numerical  integration  of  the  dynamic 
center  of  gravity  equations,  obtaining  velocity  and  displacement  components  which 
are  then  used  to  evaluate  the  equations.  After  the  evaluation  is  completed,  DYNCG 
calls  another  new  subroutine  titled  'TRANSFM'  which  updates  the  vectors  from/to 
the  new  seat/man  center  of  gravity  location.  See  Section  4. 3 for  instructions  on 
how  to  envoke  or  inhibit  the  dynamic  center  of  gravity  computation.  See  Appendix 
B for  a mathematical  presentation  of  the  dynamic  center  of  gravity  equations. 


Martin- Baker  Rail  and  Tube  Bending  Equations 

The  ICARUS  sent  ejection  model  has  been  modifled  to  simulate  either 
the  Stencil  or  the  Martin-Baker  seat  ejection  system.  The  user  of  the  model  can 
select  which  seat  is  to  be  simulated  by  setting  the  value  of  the  input  decimal  var- 
iable, 'STNCEL'  to  a rero  (0.0)  if  the  Stencel  system  is  to  be  simulated  or  to 
a decimal  one  (1.0)  if  the  Martin-Baker  system  is  to  be  simulated.  The  var- 
iable STNCEL  is  contained  in  the  decimal  portion  of  the  Default  data  file  and 
its  value  may  be  altered  by  the  user  of  the  model  via  the  normal  Card  Input 
file  procedures. 

2.3.1  Rail  and  Tube  Bending  Equations 

The  Martin-Baker  catapult  system  has  been  modelled  using  spring- 
damper  and  tube  bending  equations  to  simulate  the  forces  and  moments  acting  on 
the  seat/man  during  the  catapult  phase  of  the  trajectory’. 

Subroutine  SE  ATMAN  has  been  modified  to  interpret  the  value  of  input 
variable  STNCEL  to  ascertain  what  seat  is  being  simulated  by  the  model.  If 
STNCEL  equals  1.0  and  the  model  is  in  the  catapult  phase  of  the  trajectory  - sub- 
routine SLIPFMB  is  called.  This  new  subroutine  contains  the  Martin-B:Jier  cat- 
apult equations  of  motion.  These  equations  (the  rail  and  tube  bending)  arc 
evaluated  and  the  forces  and  moments  acting  on  the  seat/man  combination  due 
to  the  rail  interaction  with  the  slippers  and  tube  bending  are  determined  as  a 
function  of  the  velocity  and  displacements  resulting  from  the  slippers  and  tube 
bending  motions. 

If  on  the  other  hand,  STNCEL  is  input  as  a zero  or  blank  and  the  model 
is  in  the  catapult  phase  of  the  trajectory,  subroutine  SLIDFM  is  called  and  the 
forces  and  moments  attributable  to  the  Stencel  catapult  will  be  computed. 

2.4  Martin- Baker  Timing  and  Sequencing 

Subroutines  EXECUT  and  STATUS  have  been  modified  to  interpret  the 
input  decimal  value  of  variable  STNCEL  to  ascertain  which  seat  ejection  system 
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is  being  simulated  by  the  model.  As  previously  mentioned  in  Section  2.3. 1 if 
the  value  entered  into  variable  STNCEL  is  0.0,  the  Stcncel  seat  ejection  system 
will  be  simulated  by  the  model,  if  the  location  contains  the  value  of  1.0,  the 
Martin-Baker  system  will  be  simulated. 

In  the  Martin-Baker  mode  of  simulation,  subroutine  STATUS  has  been 
modified  to  check,  at  each  time  step  of  the  trajectory,  if  the  catapult  tubes  I and 
2 have  extended  beyond  the  value  specified  by  input  variable  SBRON  contained  in 
the  decimal  data  segment  of  the  Default  data  file,  if  the  tubes  have  extended  past 
this  value,  the  ‘ime  some  delta)  when  the  tube  length  equalled  the  value  spec- 
ified by  the  SBRON  variable  is  acertained.  The  values  of  the  trajectory  variables 
at  this  time  are  computed  and  displayed  and  the  seat  back  rocket  phase  is  initiated. 

Subroutine  INIT  has  been  modified  to  compute  the  time  of  drogue  pro- 
jection (TDP)  as  a function  of  the  decimal  input  variables  time  of  catapult  ignition 
(TCI)  and  the  time  required  to  elapse  (TMLAPl)  before  the  drogue  chute  can  be 
projected,  TDP  = TCI  + TMLAPl. 

Subroutines  STATUS  and  EXECUTE  have  been  modified  to  update  or 
hold  constant  at  each  time  step  of  the  trajectory  an  internal  clock.  This  clock 
has  a value  of  zero  initially  and  updating  of  the  clock  commences  at  TCI.  The 
clock  is  updated  under  the  following  conditions: 

if  the  seat/man  altitude  < ALTl 
or 

ALTl  < seat/man  altitude  ^ ALT2  and  ACCz  <TIM2GZ 
then  the  clock  is  updated  - the  integration  time  step  is  added  to  the  clock  location 
(TIMER2=TIMER2+  DE  LT  AT) . 

ALTl,  ALT2,  TIM2GZ,  and  TMLAP2  are  input  via  either  the  Default 
data  file  or  the  Card  Input  file. 

When  the  internal  simulation  clock  attains  the  value  of  TMLAP2,  the 
recovery  chute  pack  is  opened  and  the  recovery  chute  riser  line  commences 
coming  out. 

k 
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3.0 


PROGRAMMING  CHANGES 


3.1  Input/Output  Modifications 

3. 1. 1 Segmentation  of  Input  Routine 

The  original  ICARUS  subroutine  INPUT  which  reuds-In  the  aerodynamic, 
tabular,  decimal  and  integer  data  has  been  segmented  into  the  following  ten  sub- 
routines: ACT,  CONVT,  DEFAULT,  ENDCASE,  INPUT,  PTDFLT,  TABCONL, 
TABCRDS,  UPDEC  and  UPINTG.  The  control  routine,  subroutine  INPUT,  com- 
putes and  at  case  termination  outputs  the  case  execution  time,  reads-in  the  case 
title  and  user  input/output  option  cards  and  displays  them  on  the  case  output  list- 
ing. The  routine  then  calls  subroutine  DEFAULT  which  reads-in  the  decimal  and 
integer  data  from  the  Default  file.  DEFAULT  calls  subroutine  CONVT  which  con- 
verts the  decimal  data  Input  as  feet  or  radians  to  Inches  or  degrees.  DEFAUl  T 
reads-in  the  next  record  contained  on  the  Card  Input  file  and  then  transfers  pro- 
gram control  back  to  subroutine  INPLT.  Depending  on  the  values  of  QFLAG  read 
from  the  Card  Input  file,  subroutines  PTDFLT,  UPDEC,  UPINTG  and  ENDCASE 
are  called.  Subroutines  UPDEC  and  UPINTG  reads-in  from  the  Card  Input  file 
updates  to  the  decimal  and  Integer  data. 

Subroutine  PTDFLT  transfers  to  the  user's  output  file  the  updated 
decimal  and  Integer  default  data.  PTDFLT  then  calls  subroutine  TABCONL  which 
controls  the  reading-in  and  output  of  the  aerodynamic  and  tabular  data  via  sub- 
routines ACT  and  TABCRDS.  Subroutine  ACT,  reads  and  outputs  the  aero- 
dynamic data  contained  on  file  AERO.  Subroutine  TABCRDS  reads  and  outputs 
the  tabular  data  contained  on  the  Card  Input  file.  Before  transferring  program 
control  back  to  subroutine  INPIT,  PTDFLT  calls  subroutine  CONVT  which 
converts  the  decimal  data  back  to  feet  or  degrees. 

Finally,  subroutine  ENDCASE  is  called  by  INPUT  and  a check  is  made 
within  this  subroutine  to  ensure  that  all  input  data  required  to  execute  the  case 
has  been  read. 


3.1.2 


Additions  to  the  Decimal  Data 


Three  additional  input  sections  have  been  implemented  in  ICARUS, 
Sections  12,  13  and  14.  ICARUS'  input  and  output  routines  have  been  modified 
to  read  the  data  contained  in  these  sections,  from  the  Default  and  Card  Input  files 
and  to  display  this  information  along  with  the  other  decimal  data.  Section  12  con- 
tains the  input  parameters  to  the  dynamic  center  of  gravity  equations.  Section  13 
the  parameters  of  the  Martin- Baker  catapult  equations  and  Section  14  the  par- 
ameters to  the  Martin-Baker  timing  and  sequencing  logic.  For  a detailed  def- 
inition of  the  variables  contained  in  the  above  new  sections  see  Section  4. 1.2  of 
this  report. 

3.1.3  Input/Output  User  Print  Options 

An  additional  input  card  has  been  added  to  the  Card  Input  file. 

This  card  must  immediately  follow  the  case  label  card.  The  card  has  been 
added  to  the  model  input  card  stream  to  allow  the  user  of  the  model  to  select 
which  input  data  and  output  reports  are  to  be  transferred  to  the  user's  output 
file  for  subsequent  printing.  See  Section  4.1.1  of  this  report  for  a detailed  ex- 
planation of  how  this  card  is  used  to  control  the  input/output  data  display. 

3.1.4  Additional  Output  Subroutines,  Reports  and  Plot  File 

3. 1.4.1  Tabular  Displays 

At  the  request  of  Crew  Systems,  ICARUS'  output  format  was  changed 
from  a paragraph  format  to  a tabular  format  display.  A series  of  new  reports 
were  designed  displaying  the  data  in  tabular  form  (up  to  nine  components  per 
page)  with  appropriate  page  and  column  headings.  To  accomplish  this,  three  new 
subroutines  (FETCH,  NREP,  and  REPl)  have  been  implemented  and  interf.aced 
with  ICARUS'  input  and  output  routines.  Subroutine  PROUT  has  been  modified 
to  call  subroutine  REPl  when  integer  variable  IPROP,  read-in  from  the  user's 
input/output  options  card,  has  a value  other  than  zero,  indicating  that  the  user 
has  selected  at  least  one  of  the  new  output  reports  to  be  generated  and  displayed. 
Subroutine  REPl  interrogates  phase  indicator  variables  and  based  upon  the  pre- 


sent  phase  of  the  trajectory  a transfer  within  the  routine  is  made  to  the  appropriate 
Fortran  subroutine  call  statement  to  subroutine  FETCii.  The  parameters  of  the 
call  statement  indicate  the  new  reports  that  can  be  considered  for  output  in  this 
phase. 

Subroutine  FETCH  ascertains  which  of  the  new  reports  that  are  log- 
ically possible  to  output  at  this  time  have  actually  been  selected  by  the  user  to  be 
displayed.  The  values  of  the  trajectory  data  displayed  in  the  user  selected  reports 
are  then  passed,  one  report  at  a time,  to  subroutine  NREP  which  formats  and  trans-  . 
fers  the  line  of  data  to  the  associated  report  tableau  for  subsequent  output. 

3.1.4. 2 Plot  File 

An  additional  new  subroutine,  CUTTAPE,  has  been  interfaced  and  im- 
plemented in  ICARUS.  This  routine  is  called  by  subroutine  FETCH.  The  purpose 
of  subroutine  CUTTAPE  is  to  format  and  output,  to  a plot  file  at  each  step  of  the 
trajectory,  seat/man  positional,  rate  and  angular  data.  The  plot  file  generated 
by  CUTTAPE  is  then  read  by  a separate  Fortran/PLOTIO  program  via  the  Tektronix 
terminal.  The  plots  of  the  data  generated  by  this  program  are  copied  using  the 
Tektronix  hard  copying  attachment. 

3. 2 Position  Event  Time  Determination 

Subroutine  HIT  controls  the  determination  of  the  exact  time  (within  a 
given  tolerance)  in  the  trajectory  that  a specified  position  value  has  been  obtained, 
e.  g. , the  instant  in  time  when  the  catapult  tube  has  extended  its  full  length  impljing 
seat/man  separation  from  the  Martin-Baker  catapult.  /\fter  the  exact  time  of  the 
event  has  been  determined  and  output  by  a call  to  subroutine  PROUT  along  with 
the  associated  values  of  the  variables  at  that  time,  HIT  has  been  modified  to  set 
NTRY  equal  to  4 and  to  call  RKSDEO.  Entry  into  subroutine  RKSDEO  with  NTRY 
equal  to  4 forces  the  routine  to  by-pass  the  updating  of  the  work  array  which  con- 
tains the  values  of  integration  variables  at  time  ti,  ti_i,  and  ti_2.  Thus  the  values 
of  the  variables  at  time  tf  upon  entry  into  the  routine  are  not  retained.  The  rou- 
tine then  performs  the  integration  over  the  interval  from  the  event  time  back  to 
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the  time  when  the  event  was  detected  as  occurring. 

Program  control  is  returned  to  HIT  which  then  calls  PROUT  for- 
cing the  output  of  the  new  values  of  the  variables  at  the  time  of  the  event  detection. 

HIT  has  been  modified  to  update  the  Lagrange  array  YSAL  with  the  new  values  of 
the  variables  obtained  from  the  above  integration.  In  addition,  HIT  has  been  mod- 
ified (Martin-Baker  seat  only)  to  compute  the  seat  back  rocket  burnout  time  immedi- 
ately after  it  has  determined  the  time  of  seat  back  rocket  ignition. 

3.3  Simplification  of  RKSDEQ  Coding 

The  overall  coding  of  subroutine  RKSDEQ  (the  main  numerical  in- 
tegrating routine)  has  been  simplified  and  an  additional  transfer  point  was  added  to 
the  routine.  When  the  routine  is  called  with  NTRY  set  to  4,  the  updating  of  the 
work  array  before  integration  is  not  accomplished  and  the  values  of  the  current 
time  variables  are  not  retained. 

3.4  Interfacing  of  the  Martin-Baker  Time  and  Sequencing  in  Subroutine 
STATUS 

Subroutine  STATUS  has  been  modified  to  take  into  account  the  anomalies 
that  exist  between  the  Stencel  and  the  Martin-Baker  event  determination.  See  Section 

2.4  of  this  report  for  a detailed  explanation  of  the  modifications  made  to  STATUS. 

3. 5 Documentation  of  New  Subroutines 

The  following  new  subroutines  have  been  implemented  in  ICARUS; 


ACT 

ENDCASE 

RKSDEQ  (modified  and  recorded) 

CENGRAV 

FETCH 

RUNGE 

CONVT 

IDEBUG 

PTDFLT 

CUTTAPE 

INPUT 

SLIPFMB 

DANGLE 

KUTTA 

STORE 

DEFAULT 

NREP 

TABCONL 

DRICALC 

RECALL 

TABCRDS 

DYNCG 

REPl 

TRANSM 
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UPDEC 

UPINTG 

Section  4. 3 of  this  report  contains  the  detailed  flow  charts  of 
these  new  subroutines. 

3.5.1  Subroutine  ACT,  called  by  TABCONL,  reads  and  outputs  the 
aerodynamic  data  contained  on  file  AERO.  The  output  of  the  aerodynamic  co- 
efficients may  be  inhibited  by  the  user  of  the  model  via  the  iiq}ut/output  options 
card. 

3. 5. 2 Subroutine  CENGRAV,  called  by  DYNCG,  evaluates  the  equations 
of  motion  (spring  equations)  used  to  determine  the  displacement  of  the  pilot's 
center  of  gravity  from  its  static  position  during  the  seat/man  portion  of  the 
ejection  trajectory. 

3.5.3  Subroutine  CONVT  is  called  by  DEFAULT  and  PTDFLT.  Based 
tq)on  the  values  contained  in  the  locations  of  the  formal  parameters  passed,  CONVT 
will  convert  the  decimal  data  represented  In  feet  and  radians  to  inches  and  degrees 
or  vice  versa. 

3.5.4  Subroutine  CUTTAPE  generates  a file  containing  data  to  be  plotted. 
At  present  the  following  data  is  output: 

(Xs,  Ys,  Zg)  - position  of  the  seat/man  in  (EFCS) 

(Mg,  Vg,  Wg)  - velocity  of  the  seat/man  in  (SMCS) 

(((/»  9 » 9 ) “ yaw,  pitch  and  roll  angles  of  seat/man 
( 0 » 0 t ^ ) “ yaw,  pitch  and  roll  angular  rates  of  seat/man 
(Xaccel,  Yaccel,  Zaccel)  - acceleration  of  seat/man  in  (SMCS) 

3.5.5  Subroutine  DANGLE,  called  by  SEATMAN,  computes  the  derivatives 
of  the  yaw,  pitch  and  roll  angles  ( ^ , q , ,p  ). 

3.5.6  Subroutine  DEFAULT,  called  by  INPUT,  reads-in  the  decimal  and 
Integer  data  from  the  Default  file. 
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3.5.7  Subroutine  DRICALC  called  by  subroutine  DERIV  evaluates  the 
equation  of  motion  (spring  damper)  used  to  determine  the  index  of  compress!- 
bility  of  the  man's  spine  during  the  seat/man  phase  of  the  ejection  trajectory. 

3.5.8  Subroutine  DYNCG  called  by  subroutine  DERIV  computes  the 
new  location  of  the  man's  and  the  seat/man' s center  of  gravity  as  a function  of 
the  displacement  computed  in  CENGRAV  during  the  seat/man  portion  of  the 
trajectory. 

3.5.9  Subroutine  ENDCASE  called  by  INPUT  checks  to  see  that  all  in- 
put data  required  to  execute  the  case  has  been  read-in. 

3.5.10  Subroutine  FETCH  called  by  subroutine  REPl  ascertaines  which 
of  the  new  reports  that  can  be  generated  and  output  during  this  phase  have  act- 
ually been  selected  by  the  user  to  be  output.  The  user  selected  reports  are  then 
further  processed,  one  at  a time,  to  determine  which  of  the  variables  contained 
in  the  report  are  to  be  output  during  this  phase,  and  pointers  are  set  to  the  core 
locations  where  the  variables  have  been  stored  during  this  phase.  This  infor- 
mation is  then  passed  to  subroutine  NREP  which  transfers  the  line  of  data  to  the 
associated  report  output  file. 

3. 5. 11  Subroutine  IDEBUG  produces  reports  that  are  not  normally  re- 
quired to  be  displayed.  These  reports  are  used  as  aides  in  debugging  the  model 
or  to  help  in  the  analysis  of  a given  simulation  test  case. 

See  Section  4. 2. 3 of  this  report  for  instructions  as  to  how  to 
envoke  these  reports. 

3. 5. 12  Subroutine  INPUT  called  by  EXECUT  computes  and  outputs  the 
case  execution  time  at  the  termination  of  each  case,  reads  the  case  title,  user's 
Input/Output  Options  card,  and  outputs  these  two  cards  as  part  of  the  users  data 
listing.  INPUT  then  calls  subroutine  ACT,  DEFAULT,  ENDCASE,  PTDFLT, 
TABCONL,  UPDEC,  and  UPINTG.  These  routines  read  and  output  the  aero- 
dynamic, tabular,  default  and  card  input  data  files,  and  insure  that  all  data 
required  to  execute  the  problem  case  have  been  read. 
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3.5.13  Subroutine  KUTTA  called  by  DYNCG  is  the  fourth  order  Gill's 
variation  of  the  Runge-Kutta  numerical  integration  algorithm.  This  routine  is 
used  to  integrate  the  dynamic  center  of  gravity  equations. 

3.5.14  Subroutine  NREP  ealled  by  subroutines  FETCH  and  INPUT  for- 
mats and  transfers  the  line  of  data  passed  to  it  by  either  subroutine  FETCH  or 
INPUT,  to  the  output  file  designated  to  display  the  data. 

3.5.15  Subroutine  RECALL  called  by  subroutine  RKSOEQ  replaces  the 
ti  values  of  seat/man  center  of  gravity  position,  rate  data,  seat/man  vectors, 
Martin-Baker  catapult,  dynamic  response  index,  and  the  derivates  for  the  dy- 
namic center  of  gravity  and  dynamic  response  equations  with  their  corresponding 
values  at  ti_i.  These  values  are  then  used  to  determine  the  exact  time  an  event 
has  taken  place. 

3.5.16  Subroutine  REPl  is  called  by  subroutine  PROUT  when  the  new  re- 
port format  has  been  selected  by  the  user  to  display  output  data.  Subroutine  REPl 
interrogates  phase  indicator  variables  and  based  upon  the  phase  of  the  trajectory 
subroutine  FETCH  is  called  with  a formal  parameter  string  indicating  which  of 
the  new  reports  may  be  output  during  this  phase. 

3.5.17  Subroutine  RKSDEQ  called  by  subroutines  CKTME,  CKVARB,  EXECUT 
and  STATUS  has  been  modified  to  by-pass  the  updating  of  the  work  array  when  en- 
tered with  NTRY  set  to  4.  The  overall  coding  of  the  routine  has  been  simplified. 

3.5.18  Subroutine  RUNGE  called  by  subroutine  DRICALC  numerically 
Integrates  the  dynamic  reflex  index  equation.  The  numerical  technique  used  is  a 
Runge-Kutta  forth  order  method. 

3.5.19  Subroutine  PTDFLT,  called  by  subroutine  INPUT,  prints  the  default 
data.  PTDFLT  then  calls  subroutine  TABCONL  which  controls  the  reading-in  and 
outputting  of  the  aerodynamic  and  tabular  data  via  subroutines  ACT  and  TABCRDS. 

3.5.20  Subroutine  SLIPFMB,  called  by  subroutine  SE  ATMAN  during  the 
Martin-Baker  catapult  phase,  computes  the  forces  and  moments  of  the  seat/man 


due  to  the  rail  interaction  with  the  seat  sliiq>ers  and  catapult  tube  bending. 


3.5.21  Subroutine  STORE  called  by  subroutine  INIT  and  RKSDEQ  saves 
the  values  at  time  ti  of  the:  seat/man  center  of  gravity  position  and  rate  data, 
Martin-Baker  catapult  data,  DRI  data,  and  the  deriviatives  of  the  dynamic  cen- 
ter of  gravity  and  reflex  index,  after  each  integration. 

3.5.22  Subroutine  TABCONL  called  by  subroutine  PTDFLT  controls  the 
reading  of  the  aerodynamic  and  tabular  data  via  subroutines  ACT  and  TABCRDS. 

3.5.23  Subroutine  TABCARDS  called  by  subroutine  TABCONL  reads  the 
tabular  data  from  the  Card  Input  file  and  transfers  the  data  to  the  output  file. 

3.5.24  Subroutine  TRANSFM  called  by  subroutine  DYNCG  updates  the 
vectors  from/to  the  seat/man  center  of  gravity  during  the  ejection  phase  of  the 
trajectory. 

3.5.25  Subroutine  UPDEC  called  by  subroutine  INPUT  reads  from  the 
Card  Input  file  updates  to  the  decimal  data. 

3. 5. 26  Subroutine  UPINTG  called  by  subroutine  INPUT  reads  from  the 
Card  Input  file  updates  to  the  integer  data. 

Section  VI- 3 of  this  report  contains  detailed  flow  diagrams  of 
the  above  new  routines. 

4.0  USER  INPUT  DATA  INSTRUCTIONS 

4.1  Definition  of  Input  Variables 

4. 1. 1 Definition  of  User's  Input/Output  Options  Card  Format  (612,  3X, 
3511). 

Subroutine  INPUT  has  been  modified  to  read-in  an  additional  in- 
put card.  This  card  is  mandatory  and  must  immediately  follow  the  case  title 
card.  The  card  contains  the  following  user  coded  output  display  switch  settings; 
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Card 

Col. 

2 


4 


6 


8 


12 


Col. 

Name 

IPROP  A zero  or  blank  entry  will  force  the  program  to  dis- 
play output  variables  In  original  ICARUS  format. 

A '1'  keypunched  In  this  column  will  force  the  pro- 
gram to  display  output  variables  In  the  new  report 
format. 

Reports  2-23 

A *2'  keypunched  In  this  coU:mn  will  force  the  pro- 
gram to  display  output  variables  In  both  the  original 
and  new  formats. 

IPLOT  A zero  or  blank  entry  In  column  4 of  the  card  will  In- 
hibit the  generation  of  the  file. 

A '1'  keypunched  in  column  4 of  the  card  will  force 
the  program  to  generate  a file  containing  data  to  be 
plotted. 

IRECH  A zero  or  blank  entry  in  column  6 of  the  card  will  In- 
hibit the  outputting  of  the  recovery  chute  angle  and 
reference  area. 

A keypunched  in  column  6 of  the  card  will  force 
the  program  to  display  the  recovery  chute  angle  of 
attack  and  reference  area  after  each  iteration. 

lAERO  A zero  or  blank  entry  in  column  8 of  the  card  will 
inhibit  the  outputting  of  the  seat/man  aerodynamic 
data. 

A ’1'  keypunched  in  column  8 of  the  card  will  force 
the  program  to  output  the  aerod3mamic  data. 

IDFLT  A zero  or  blank  entry  In  column  10  of  the  card  will 
Inhibit  the  outputting  of  the  default  data. 

A *1’  keypunched  in  column  10  of  the  card  will  force 
the  program  to  output  the  default  data. 

ITABP  A zero  or  blank  entry  in  column  12  of  the  card  will 
Inhibit  the  outputting  of  the  tabular  data. 
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A *1*  keypunched  in  column  12  of  the  card  will  force 
the  program  to  output  the  tabular  data. 

Card  Col. 

Cols.  Name 

16, 17, 18, ...  35  IRSW  (1, 2. . . 23)  A 'O'  or  blank  entry  in  these  columns  will 

force  the  program  to  inhibit  the  printing  of 
corresponding  new  reports  (Reports:  1,2,3 
23). 

A '!'  keypunched  in  these  columns  will  force 
the  program  to  print  the  corresponding  new 
reports  (Reports:  1,2,3. ..23). 

Presently,  reports  2 to  23  have  been  defined  and  implemented. 
Subroutines  ICARUS,  INPUT,  PROUT  and  RECOV  have  been  modified  and  three 
new  subroutines  (FETCH,  NREP,  REPl)  have  been  added  to  the  model  to  inter- 
pret the  above  inputted  switch  settings  and  to  generate  the  user  requested  output 
reports. 

4.1.2  Definition  of  the  New  Decimal  Data  Sections 

The  following  three  additional  input  sections  have  been  added  to 
the  decimal  data  segment  of  the  Default  file: 

4. 1.2.1  Section  12,  Dynamic  Center  of  Gravity 

Section  12  contains  the  following  12  input  parameters  to  the  dy- 
namic center  of  gravity  equations: 

CX  XSLACK  SXP 

SXN  CY  SY 

CZ  ZSLACK  SZP 

ZBOT  SZNl  SZN2 

Where:  CX,  CY,  and  CZ  are  the  damping  constants  in  the  X,  Y and 

Z directions  respectively,  expressed  in 
units  of  pcr-second. 

SXP,  SXN  spring  modulus  constants  in  the  X direction, 

expressed  in  units  of  per-second 
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XSLACK 

dead  zone  in  X direction  (0<  X S XSLACK) 
where  the  frictional  acceleration  term  (FX) 
is  set  to  zero. 

SY 

spring  modulus  constant  in  the  Y direction, 
expressed  in  units  of  per  second  2. 

SZP,  SZNl, 
and  SZN2 

spring  modulus  constants  in  the  Z direction, 
expressed  in  units  of  per  second  2. 

ZBOT 

bottoming  distance  in  the  Z direction,  ex- 
pressed in  feet. 

ZSLACK 

upward  slack  in  the  Z direction  (0<  Z < 
ZSLACK)  where  the  frictional  acceleration 
term  (FZ)  is  set  to  zero,  expressed  in  feet. 

Reference  Appendix 

B for  a mathematical  presentation  of  the 

above. 

4.1.2. 2 Section  13,  Martin-Baker  Catapult 

Section  13  contains  the  following  26  input  parameters  to  the  rail 
and  tube  bending  equations  plus  the  program  control  switch  STNCEL  indicating 
which  ejection  sect  system  is  to  be  simulated  by  the  model  (STNCEL  set  to  0.0 
implies  the  Stencel  system,  set  tc  1.0  implies  the  Martin-Baker  system): 


XKSLIP 

ZKSLIP 

XCSLIP 

ZCSLIP 

UMSLIP 

XYKTOR 

XYCTOR 

YKTUB 

PKTUB 

RKTUB 

YCTUB 

PCTUB 

RCTUB 

UMTUB 

TUBLTl 

TUBLT2 

TUBLT3 

ELMOD 

DIAMIl 

DIAMOl 

DIAMI2 

DIAM02 

XCTCP 

YCTCP 

ZCTCP 

STNCEL 

RAILNT 

IMiere: 

XKSLIP,  ZKSLIP 

are  the  spring  modulus  constants  in  the  X 
and  Z directions  of  the  rall/slipper  equations 
in  units  of  /sec 
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XCSLIP,  ZCSLIP  are  the  damping  constants  in  the  X and  Z 
directions  of  the  rail/slipper  equations  in 
units  of  /sec. 

UMSLIP  coefficient  of  friction  used  in  rall/slipper 

equations  (Y  direction). 

XYKTOR  rail/slippcr  torsional  spring  modulus  con- 

stant in  the  pitch  plane  PC-Y)  in  units  of 
(ft-lbs/rad). 

XYCTOR  rail/slipper  torsional  spring-damper  con- 

stant in  the  pitch  plane  pc-Y)  in  units  of 
(ft-lbs  - sec/rad). 

YKTUB,  PKTUB,  catapult  torsional  spring  modulus  constants 

RKTUB  in  the  yaw,  pitch  and  roll  direction  in  units 

of  (ft-lbs/rad). 

YCTUB,  PCTUB,  catapult  torsional  spring  damper  constants 

RCTUB  in  the  yaw,  pitch  and  roll  directions  in  units 

of  (ft-lbs-sec/rad). 

UMTUB  catapult  coefficient  of  friction  in  Y direction. 

TUBLTl,  TUBLT2,  length  of  tubes  1,  2 and  3 respectively  in 

TUBLT3  feet. 

DIAMIl,  DIAMI2  inside  diameters  of  tubes  1 and  2 respec- 
tively in  feet. 

DIAMOl,  DIAM02  outside  diameters  of  tubes  1 and  2 respect- 
ively in  feet. 

XCrCP,  YCTCP,  coordinates  of  the  seat  contact  point  with 

ZCTCP  catapult  tube  in  feet  (CCS). 

RAILNT  length  of  slipper  rails  in  feet. 

See  Appendix  C for  a mathematical  presentation  of  the  use  of 
the  above  variables. 
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4. 1.2.3 


Section  14,  Martin-Baker  Phasing 


Section  14  contains  the  following  7 input  parameters  to  the  Martin- 
Baker  drogue  chute,  seat  back  rocket  ignition,  recovery  chute  and  seat/man  sep- 
aration simulation: 


TMLAPl 

ALTl 

SBRON 

Where:  TMLAPl 

TMLAP2 


ALTl 


ALT2 

TIM2GZ 


TMLAP2  ALT2 

TIM2GZ  SMSEPF 


is  the  time  (from  catapult  ignition)  of  drogue 
projection  for  the  Martin-Baker  system. 

is  the  time  that  the  recovery  chute  clock  must 
advance  to,  from  a zero  initial  value,  be- 
fore the  Martin-Baker  recovery  chute  pack 
will  open. 

altitude  below  whieh  the  recovery  chute  clock 
is  updated  automatieally  until  T..ILAP2  is 
reached  (in  feet). 

if  seat/man  are  above  this  altitude,  the  re- 
eovery  ehute  eloek  is  not  updated. 

vertical  acceleration  (in  G's)  which  will 
stop  the  clock  from  being  updated  if  the 
man's  altitude  is  below  ALT2  and  above 
ALTl  and  his  vertical  acceleration  is  great- 
er than  TIM2GZ. 


SMSEPF  resultant  force  that  the  recovery  chute  must 

exert  (greater  than  or  equal  to)  on  seat/man 
to  cause  the  man  to  separate  from  the  seat, 
in  units  of  lbs. 

SBRON  distance  seat  must  travel  up  the  catapult  tube 

before  the  seat  back  rocket  will  ignite,  in  feet. 

Changes  to  the  decimal  data,  defined  above  for  Sections  12,  13  and 
14,  can  be  accomplished  at  run-time  via  the  normal  Card  Input  procedures,  i.e., 
entering  in  the  decimal  data  section  of  the  file:  the  variable  name(s)  (columns  2-7) 


I 

I 


I 

I 

i 

f 

i 


;! 


-20- 


28-33  and  54-59),  and  corresponding  values  (columns  8-27,  34-53,  and  60-79) 


Reference  Section  3 of  NWL  Technical  Note  TN-K  7/74  dated 

February  1974. 

4.1.3  Definitions  of  the  New  Integer  Data  Variables 

The  second  segment  of  the  Default  data  file  contains  the  integer 
input  data  to  the  model.  The  first  word  of  this  segment  contains  the  count  (30) 
of  the  number  of  integer  data  values  that  immediately  follow.  The  next  two  re- 
cords contain  the  integer  values  to  be  read-in  under  a 1615  format.  The  first 
four  words  of  the  30  word  array  were  used  by  the  original  implementors  of  the 
model.  The  remaining  26  locations  were  set  aside  for  future  use. 


In  the  revised  version  of  ICARUS  5 of  the  26  extra  words,  words 
5 through  9 of  the  30  word  integer  array  INT,  have  been  used  as  follows: 


Word  Card  Columns 


6 

7 

8 


21-25  Set  by  user  of  the  model  to  control  the 

computation  of  the  dynamic  center  of  grav- 
ity. 

0 - entered  in  this  word  (via  the  Default 

file  only)  will  force  the  model  to  by- 
pass the  dynamic  center  of  gravity 
computation. 

1 - entered  in  this  word  will  force  the  model 

to  compute  the  dynamic  center  of  grav- 
ity. 

26-30  Set  by  the  user  of  the  model  to  request  add- 

31-35  itional  output  reports  not  normally  displayed  - 

36-40  model  debugging  or  case  analysis  aides. 

Entering,  via  the  Default  file  only,  the  follow- 
ing values  into  these  words  will  force  the 
model  to  generate  the  following  reports: 

1 - Prints  all  forces  and  moments  on  seat/ 

man. 

2 - Prints  values  cont:iined  in  locations  NTRY, 

NUMEQS,  DELTAT,  DELTARP,  WORK 
(1-5). 
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3 - Prints  forces  and  moments  due  to  catapult. 

4 - Prints  Martin-Baker  catapult  displacements 

and  rates. 

5 - Prints  the  Quaternions. 

6 - Prints  CTLWDS  common  block. 

^ 41-45  Set  by  the  user  of  the  model  to  control  the  gen- 

eration of  the  Tektronix  plot  file. 

0 - entered  in  this  word  will  inhibit  the  gen- 

eration of  this  file. 

1 - entered  will  force  the  model  to  generate 

the  file. 

4. 2 Program  Load,  Execute,  Output  and  Load  Map  Procedures 

4.2.1  CALL,  PRLD  - this  procedure  loads  the  ICARUS  program  and 

associated  data  files  and  executes  the  program. 

KRONOS  PRCX:EDURE  PRLD: 


00110  ir»N.MS«ouTPiiT. 

00111  Rrt'lND«OUTPIIT  .TAPP  1 1 • T ftPE  1?.  T APp;  1 T APF  1 A . T APt  1 5 . T APp  1 ^ , T A^E  1 7 . 

00112  REPlNn,TAPElR.TAPp  If^.T/iPEPn.TAOKZl  , T APE . T ApF?3 , T APp ?a  . T Aap 25 . 

00113  9E'-»IN0*TAPEPf'.TAPP?7.T  Ade2R.TAPE?P»  TAPr''n.TAPE31  . T APE  32  • T APE.  33  . 

001 1 A R^WI^JD•TAPE3A.TAPE^S*TADE3ft,TAPE37»TAPE3B,TAPE39»TAPEA0.TAOEAl , 
00115  RF'*»IajD.TAPP4?.TAPFA3.TapEAA.TAPEAS. 

00120  GET»TAPFl=DPnATA. 

00130  GFT,CARnT\. 

OOIAO  GET»TAPF2=AER0. 

00150  GET. MEW. 

00155  GET*COMrLIR/l)M  = SYSTFM, 

00160  GFT»maTl  IR/iJM  = SySTc-m. 

00163  Rn  *127777, 

00170  LP^ET.L rP=MATLIR. 

00171  Ln«:ET,LTR  = COMCLIH. 

OOIBO  LOAD. NEW. 

00190  EXECUTE.  ICAPIIS.CARhin. 

00200  EXIT. 


4.2.2  CALI,  TAPOUT  - this  procedure  copies  the  new  report  files  onto 

the  output  file. 
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A 


i 

j 

i 


1 


s 

i 


KRONOS  PROCEDURE  TAPOIH*; 


REV/lN0»TAPEll,TAPE12,TAPf:i3»TAPE14,  TAP£15.TAPE16.TAP£17,  r6PEl3»TAP£19. 
.REWIN0»TAPE21  -.TAPE22»TAPE23»TAP£2AtTAPEi?t>»TAP£2t>,  rAPE27»TAPc23»TAP£29, 
REWIND*  TAPE20* TAPE30»TapE31 , rAPE32* TAPE 33. 

REWlND»TAPE3C.-TAPE35-*TAPE36*TAPE37»TAPE3d*TAPE39,TAPEA0*TAPEAl .. 

REWIMD»TAP£42*  TAFEA3*TAPE^4,TAPEA5. 

COPYBF*TAPEl 1 . 

C0PYBF*TARE12. 

C0PYBF»TaPE13. 

COPYBFtTAPElA. 

-COPYBF*  tape  15. 

C0PYBF*TAPE16. 

C0PYBF»TAPE17. 

- C0P-YBF»TAP-£33 ^ 

COPYBF, TAPE30. 

COPYBF, TAPE31. 

-COPYBEL,  TAPE3^ 

COPYBF, TAPE18. 

COPYBF, TAPE19. 

_COP-YBF,TAPE20. , 

COPYBF, TAPE21 . 

COPYBF, TAPE22. 

...COPYBF.,  TAPE23. 

COPYBF, TAPE29. 

COPYBF, TAPE25. 

-COPYBF,  TAPE26-. : 

COPYBF, TAPE27. 

COPYBF, TAPE28. 

-COPYBF,  I APE29..- : 

COPYBF, TAPE3A. 

COPYBF, TAPE35. 

..C0P.YBE.,TAP£36_ 

COPYBF, TAPE37. 

COPYBF, TAPE38. 

..COPYBF,  TAPE39.— 

COPYBF, TAPEAO. 

COPYBF, TAPEAl. 

COPYBF,  TAPE.42. 

COPYBF, rAP£43. 

COPYBF, TAPEAA. 

.C0P.YBF,TaP£A5. 


L 
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4.2.3  CALL,  MAP  - this  procedure  loads  the  ICARUS  program  and 

associated  data  files  imd  generates  a full  load  and  subroutine  cross  reference 
map  of  the  program. 


KRONOS  PROCEDURE  MAP; 


get»tapei=dfoata. 

GET  »C0MCLIB/UN  = SYST£M. 
GETtCAROlti. 

GET*-TAPE2=AERa. 

GET* NEW. 

GET  »MATlIB/UN=SYSTEM. 

RFL»  1277-77 

ASSlGN»MS.'*jDj£rPUT.  

LD5ET»Li8=MATL l«»MAP=SdEX . 
-LDSET»L13=C0MCLIB. 

load»mew»matlib. 

NOGO. 

DISPOSE.OUTPUT=PR/F1=CD3002.. 


5.0  HOW  TO  INSERT  NEW  DATA  AND  NEW  REPORTS 

5. 1 Additions  to  the  Default  Decimal  Data  (Segment  One  of  the  File) 

Whenever  additional  decimal  input  parameters  are  required  by  the 
model,  i.  e. , extensions  to  the  model  such  as  the  dynamic  center  of  gravity  equa- 
tions, a new  data  section  can  be  appended  to  the  last  section  (currently  Section 
14)  by: 

• Changing  the  first  record  (one  word  in  length,  15  format)  of  the 
Default  data  file  (DFDATA)  to  reflect  the  number  of  additional 
decimal  words  to  be  input.  Currently  this  word  contains  the 
Integer  value  of  253.  If  a new  section  were  to  be  added  con- 


L. 
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talning  six  additional  variables,  the  word  count  would  be  in- 
creased to  259.  The  values  of  the  six  new  variables  to  be  read- 
in  would  be  inserted  immediately  after  the  253rd  word  (words 
254,  255, . . . 259)  in  a 3E20. 14  format. 

• Subroutines  INPUT,  INIT,  PTDFLT,  UPDEC,  DEFAULT  and 
the  appropriate  labelled  common  blocks,  dimension  and  integer 
array  must  be  updated  to  reflect  the  additional  variables. 

The  253  (current)  decimal  data  values  are  displayed  on  the  user's 
output  file  in  the  14  data  sections  as  follows: 


Section 

Begin  Word 

End  Word 

Length 

1 

1 

12 

12 

Angular  data 

2 

13 

85 

73 

Positional  data 

3 

86 

97 

12 

A/C  data 

4 

98 

103 

6 

(Seat  back  rocket  angles) 

5 

104 

108 

5 

Slider  block 

6 

109 

111 

3 

Dart 

7 

112 

145 

34 

S/M  inertia,  weights,  area 

8 

146 

163 

18 

Chute  data 

9 

164 

174 

11 

Chute  data 

10 

175 

183 

9 

Event  time  date 

11 

184 

207 

24 

Event  time  steps 

12 

208 

219 

12 

Dynamic  C.G. 

13 

220 

246 

27 

M.  B.  catapult 

14 

247 

253 

7 

M.  B.  phasing 

5. 2 Additions  to  the  Default  Integer  Data 

(Segment  2 of  the  File) 

The  first  record  after  the  decimal  data,  1 word  in  lenfrth  (15-format) 
indicates  the  number  of  integer  words  that  immediately  follow.  Currently  this 
count  is  set  at  30.  The  next  two  records  contidn  the  integer  data  values  in  a 16 
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15  format.  The  first  4 words  of  the  30  word  integer  array  have  been  formally 
defined  by  the  original  implementors  of  the  model  and  have  been  equivalenced  to 
the  following  integer  data  names:  ISMA,  ISTOP,  ISMS  and  ISEAT  respectively, 
see  pp.  - 63,  71,  72  of  the  Program  Maintenance  Manual  (Reference  2)  for  a de- 
tailed definition  of  their  use  in  the  model. 

The  first  4 integer  data  inputs  may  be  updated  via  the  Card  Input 
file.  The  remaining  26  words  have  not  been  formally  defined  for  updating  using 
the  Card  Input  file  procedure.  The  26  words  can  currently  only  be  modified  by 
changing  their  values  in  the  Default  file.  Words  5 to  9 are  currently  being  used 
for  plotting,  program  case  analysis  and  debug  purposes. 

See  Section  4. 1. 3 of  this  report  for  further  details  regarding  their 
use  in  the  model. 

5.3  Implementation  of  New  Reports 

The  first  card  immediately  after  the  case  title  card  is  the  user  INPUT/ 
OUTPUT  Display  option  card.  This  card  contains  35  switches  (one  for  each  report) 
of  which  23  are  currently  in  use.  These  switches  are  set  by  the  user  to  control 
the  printing  or  non-printing  of  each  of  the  35  reports.  If  a new  report  is  to  be  added, 
switches  24  to  35  are  available  to  be  assigned  to  the  report.  Report  page  and  column 
headings  can  be  entered  by  inserting  the  appropriate  page  and  column  headings  in 
the  data  statement  array  allocated  in  subroutine  NREP.  Subroutines  RE  PI  and 
FETCH  must  be  updated  to  call  the  new  reports  with  formal  parameters  pointing 
to  the  data  locations  to  be  printed  and  the  format  to  be  used  in  displaying  the  loc- 
ations. 
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6.0  ICARUS  FLOW  DIAGRAMS 

6. 1 Macro  Structure  of  ICARUS 

Control  Section 

Overall  model  data  initialization. 

Case  study  data  initialization. 

A/C  and  seat/man  requests  for  trajectory  computations  and 
output. 

Adjustments  to  the  trajectory  time  steps. 

Termination  of  case  studies. 


Event  Control 

Monitors  the  status  of  position  and  time  events. 

Initiates  output  request  for  message  indicating  a particular 
event  has  occured. 

Initiates  requests  for  the  determination  of  the  exact  time  that  an 
event  has  occurred  and  to  output  this  time  and  the  values  of  the 
trajectory  variables  at  that  time. 

Dynamic  scheduling  and  initialization  of  the  equations  of  motion 
to  be  evaluated  at  each  phase  of  the  trajectory. 

Controls  overall  time  step  ( A t)  to  be  used  during  each  phase. 


;;ntir.  UKSlII'.f} 


iubr.  IJNEAH  | f lUioil  I I Subr.  SKNT 


-3 


COnPUTE  OER[VflTlVES  OF  YRH.  PITCH.  RND  ROLL 


-6.2631653071795  ♦ PHI 


£N0 


6 


CONI.  ON  PO  3 

-48- 


po  2 Ql 


r 


( 


DYt;cc~| 


1 C9LL  CtNGf.ftY  1 

1 

1 

12 

1> 

7 

1 COMPUTE  RCGM  | 

XCG:1=XCGM0  + DXCC-M 
YCC-M=YCGM3  t DYCC-M 
ZCGMriCGMQ  •(■  DICGM 

12 

7 

1 COMPUTE  CCCrGM  | 

7 

XCGSIlUCiiMsXCGM  t MSRcXCGS)/MS 
YCG3l-il  = C;incYCr-M  + ilSPcYCGGl/tlS 
ZCGGMl  = ( i'r:i'::ZCG:i  t MSfirZCGS  )/HS 

t 

12 

[cnn^ui 

, 

E OPCCf  1 

y 

DXCC-=XCGSMl  - XCGStI 
DYCGrYCCCin  - YCGSM 


Dicrrr.zca'itu  - rcGzn 

FORnmi  lOX  .££14  .G.SHx-rc  ) 
FORt'.mc  lOX.GCU  .S.SX.SHiiCCSM./ 


CONT.  ON 
~r>\- 


nfCCSM= 

1 ZCDSM- 

rCGSMl 

jrCGGttl 

1 

CRLL  TRnKVlKQY)  | 


RETURN 


END 


eiLS_rir)rL_ 


or_ 


-r.s- 


ec_r'_tLt’nL_ 


SIJSnOUTlK'S  rfTCHCK'f.'SPOR'T.Ul  .IJ2.!n.Ky-.W5.K’S.!r7.N8.N9.NI0 
1 Hll  .W12  .t)l3.in4,l>!lt3.IKG.(n7.HlG  .!.U9,W20) 


COfM.  ON  ro  2 
-;■>!>- 


EO_l cj: IQ 


CONT.  GM  PO  n 
-(ifi- 


ao-iO-C 


-7(1- 


COHT.  Or4  PO  13 


or*  to  nr 


CRLL  NREr-ClQ.? . WORK . [ S I , FtRCS . 1 0 .Ki-iOEHT . 10 ) I 


CONT.  ON  ro  14 
-71- 


PC.J1.CF 


com.  on  po  15 

72-  E{LJLl_nf 


IS2=.34->-IUPC  IS2=»3I-i-IUP  ] 


CRLL  »P'EP(  15.6.‘^»PH.IGl.»aKK.IS2) 


ICO  TO 

300 

r ISlaZlflUpS 

IS2.:rT3fLUP 

.A 

7 

1 CP.LL  NP£P(  l‘i.S.t^OR:t.IOl.’nGRK.IS2) 

jJ 

, \ 

7 

|G3  ro 

300  > 

-7.3- 


CO»T.  ON  PO  16 


P.G..l5_0r 


I G3  TO 


CALL  NREPt  l6.6.r3RC3.2Z.R«0£IJr.22J 

f G3  TO 

7 

"03  > 

COLL  NREPC17,S.Fe 

RCS.25.Rtt0EKr.25) 

^ ^ 

03  TO  sea  ^ 


ZOCC  = I.RCC/5PG.OSO 

__  _ . ■ A 

7 

I COLL  N‘?E?(  rj.G.n.RL  .1 ) 

7 

1 03  TTJ  SC3-3> 

CTiLL  NREPl  19  .*7  .DXCGM . I .XHCCO. . I .TXFRCE . 1 1 


1 


V 

COHT.  0 4 ro  10 


PCL_l.7._Ci; 10 


-7G- 


efi_iaj3E 


ST2( I )=.5857aB643762690rGTl( [ U . l2l32Q3435S964eST2( I ) 1 

[i] 



1 XC  n=X[  I ) + l .70710G781  185Sk(ST1  { n-5T2(  n ) | 

X 

7 

1 ST2(  I J=3.4U2l3562373l«STl(  [ J-4 . 12l3203435596cSr2(  I J | 

S 



7 

5 > 

< \ 


I END 


-7<)- 


ea_2_EiNi.L_ 


CONT.  ON  PO  2 
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ec.  l__QF. 


3 


e£L2._QF__5;_ 


39 


URITE(  IRNl  .003)  CPR0liDS(  I ) .1  = 1 .8)  .IPfiGC  IKN)  | 

"■■■  — ■ . - - - ■ . ■ J 


HRlYErrRNr.Q04)  HHEN.IRN 

FORURTCIH  .SOX.SHDflTEt.HlO.^X.THREPORTC .I2.3H):  ) 


, 

RRITEC IKN2 .JFMTC 1 .IRN))] 


’^7 

GO  TO  12  '> 


IF(  ICOUNK  IRt4)-52)  12.12.11 


IPRGC IRN)= 

IPRGC IRN)+1 

■ J 

7 

WRITEC IRNl  .803)  (PROEDSC 1 ) .1  = 1 .8) .IPflOC IKN) 
FORRflrC  IHl  .IH  .8fll0.30X,KPflGEt!.  13.//) 

J 

7 

ICOUNK  IKN)=0 

K'RITEI  IRKl  .804)  WHEN.  IRN 
WRITEC IKNl .JFMTC 1 .IKN) ) 

ICOUNTC IRN  ) = 

ICOUNTC IRN)+1 

CONT.  ON  PO  5 
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E£}_4__C£. S 


r 

I 

I 

39 

jltRirEcTRi^r.NFfifT  r?fiE  .Tpc  r j . i=i  ,nvTT 


«RrTE(  IRin  .IFMK  1 .K3)  TI«E.(P(  I J .1  = 1 ,I-!V1 ) 


I 


PRIKT  5« 


IF( .N0T.(J.£Q.7  .OR.  J -EQ.  13)) 


GO  TO  100 


QPGCRT^^PGCKT^l 


PRINT  6Q.CPGCNT 


PRINT  54- 
PRINT  64-.CSXCNT 


PRINT  55.  (80tNR^iC  I ) .GLCCKK  I ) .laT.NBl) 


loo  TO  90 


PRINT  55.  (BOP-NRK(  I ) .BL0CK2(  I ) .I=a.KB2) 


GO  TO  90 


CONT.  ON  PO  3 

-8()- 


eO.  2.  -Qg 


PRINT  SS.  (803NRt;(  n. BLOCKS ( n .r=:l.NB3)  | 
TgB  TQ^  90  ^ 


PRINT  55.  (B04Nfif;(  r).SL0CH4(  1)  | 

I 2 ^ 

I 03  TO  9Q 


PRINT  55.  (BGSKRt^I  n.BLCCKSC  n .I=:I.NB5)  | 

t V- 

|GC  TO  90  ^ 


PRINT  55.  (BGSNRMID.BLOCKSn)  .t^l.NBSTl 

‘ " —I 


foo  TO  90  t> 


PRINT  55.  (607NRi1(  n.eL0CK7m  .I=l.KS7)  | 


00  TO  90  t> 


CONT.  ON  PO  4 
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1 PRINT  S5.(Bl3Nfif1(  n.aL0CKl3(  I).I  = l.N3t3y| 

I go  to' 

SO  > 

1 PRINT  55.(Bl+URM(  n.BLOCiTUI  n.I:a,NP>l4)  | 

, 

, 2___ 

I GO  TO  £5  ^ 


osxccr  =>osxcNr  i 

QPGCMT  =^PPGCNr-^l 

\ 

7 

PRINT  sa 
PRINT  5G 

I.8PGCNT 

39  i 

7 

IPRINF  THE  CONTROL  MOROS  | 

« 5 i 

^ T 

1 PRINT  57.(INrNflh(n.INT  (n.I=l.NtN)  | 

^ 

1 ■wtrFP  PV1RRM£TERS~^D  CONTROL  KQKDS  RR£  PRCNTfO****  | 


^ 

CONVERT  PLL  PflKRhETERs'cQNrOtNtNG  PHGLES  IN  D£0k~££S  TO  1 


if 


CONT.  ON  PO  n 
-80- 


PP-5.  or  6 


eQL_£_JLINCL_ 


I 


CONTINUE 


xsrcp  = SP.VETCPC  INITIRL.  I ) 
YSTCP  = GRVETCf'C  INITIRL. 2) 
ZSTCP  = SflVETCPI INITIAL. 3) 


JET 


I?. _ 

VALUES  PON  SLIPPMfi 


I 


CONT  . ns  , 


< 00  400 

1=15.22  > 

! 

t 

RUG«(I)  = SRVESLPI INITIAL. 1-14) 
CONTINUE 

f 

< 00  500 

1=34.45  > 

5Z 

RUGKdJ  = 8flVESLP(rNITIRL.I-33«-8) 
CONTINUE 


, 12_  , 

r^ESET  VRLUES  FOR  ORICRLC  | 


< 00  600  I=ll 


V 

RUGKCI)  SRVEORK  INITIRL.I-IO) 
CONTINUE 


L2__i , 

RESET  OERIVRTIVES  FOR  OYKCQ  RNO  ORICRL^ 


< DO  70Q  _I=82.89  > 


^ 

RUGKd)  = SfiVEOYI  INITIRL.I-81 ) 
CONTINUE 


CONT.  ON  PQ  4 
-93- 


eiLJ_0£_ 


WRITE( 11 .1000)  1SMITCH.8RVECC( IHIT IRL . 1 ) .INIT IRL  . IFINRL 
FORMOTC  IHO  .5X  .SHRECflLL  .5X  .8H1SM 1 TCH= , 1 3 .3X  .SHXCGStlr  .El4  .5. 


i* 


-‘)1- 


PP  4 FIN’PL, 


FCSTNCEL  .EQ.  0.  .OR.  ISLIO  .EQ.  1) 


r 


<rb6  730  .O 


Qo  730  J=l  .T~> 


I Fftfinnc j. £ ) =FH(£.j)  I 


rCflLL  FETCH  £ 4 .20,21. 22.23J  | 

[iizzz: 


:£ 


[ go  to  ~74  ^ 


CONT.  ON  PO  3 
-I)G- 


eo  ^._0F. 


CONT.  ON  PO  6 
-1)5)- 


e0_5..0P Q 


IF(  IHORD-l  J12S.U7.U7 
IFdSC-n  U8.U9.12S  I 

C9LL  FETCH(2.U.12J  ' 

IFCIRECOV  -2J  1181.1182.1182  • 
CRLL  FETCH!  l.UJt  GO  TO  12S 
CALL  FETCH! 1. IS) 

IF!  I6f1SEP)120.120.12l 
CRLL  FETCH! 1.161 


eQ..7__QE_  J 


-102- 


ep  a-_eiN9L_ 


RKSOEQ 


4 

SUBROUTINE  RKSOECIQOERIV) 
IPHflSE  = NTRY  » 1 


4. 


m NTRY 0 

CONTINUE 

rNPPSI 

Z , 

, 39  , 

riNiTIRLlZRflON  ROUTINE  | 


OYN(  I )=S' 

IQ( 1 I=QYN( 1 1 


, 39 4. , 

I POINTERS  IN  WORK  RRER  | 


CONT.  ON  PO  2 
-103- 


eO  I .OF 


QYN(2«n=ia(ui) 

OYPNl  r J=(l4-9)*NUflEQS^S 

B 

CJ 

PI 

1 

? 

aKN=IQ(4) 

0 

fi2 


OKENO  s QQPN  «•  NUMECS  - I 


OI=GYN(  D^NUnEQS 


00  15  QJ  =:  QI.QKENO 


M0RK(QJ)=0. 


COMPUrE  INiriRL  VALUES  OF  DERIVRriVES 


CALL  0bERlV(A0~RK'(2 ) .NUtiEOS  .AORKC 5 1 .MORKC  ODYN  m 


NTRY  = I 


RETURN 


CONT.  ON  PO  3 
-101- 


I 


eo.2  OF . - 11 


P0.3__0F___l 


CONT.  ON  PO  S 
-lOG- 


ee-  +__0F_  _u__ 


efi_  5 -jf 1 1 


CONT.  ON  PO  7 
-108- 


Pa^__DE Ll_ 


I CONFINUE 


CONT.  ON  PO  8 
-109- 


EO  7__P£ 


foo  TO 

1 

^20^ 

1 H0RK(LJ=H0RK(4.)*l  .707106781  186S»!(W0RK(J)-H0RK(N))  | 

_ 1 

_ 

1 «0RK(N)=3.4142l3562373l«K0RKCJ)-4.l2l32a3+35S96*K0RK(N)  | 

\ 

7 

I GO  TO  320~^ 


CONT.  ON  PO  10 
-111- 


PO  9 


-113- 


eoji. Finnic 


RUNGE 


INTEGER  FUNCTION  RUNGE  (N.Y.F.X.H) 


39 i 

THIS  SUBROUTINE  IS  COPIED  VERBATIM  FROM  PAGE  374  OF  CRRNflHflN , 


LUTHER,  AND  WILKES:  APPLIED  NUMERICAL  METHODS.  WILEY.  1969. 

■ EKKKKKKMKatKUKKKKKKKaKKKKKHKKHKMKXXKKKKKKKKHKaKaiKKKKKKKKaKKKKKaKIIK 


THE  FUNCTION  RUNGE  EMPLOYS  THE  FOURTH  ORDER  RUNGE  KUTTA  METHOD 




[ 

WITH  RUNGE"S  COEFFICIENTS  TO  TO  ! 

INTEGRATE  A SYSTEM  OF  N SIMULTAN-  I 

EOUS  FIRST  ORDER  ORDINARY  DIFFERENTIAL  EQUATIONS  F( JJ=:OY( J)/DX , 


Rjrl .2 . . . . .N ) . ACROSS  ONE  STEP  OF  LENGTH' H IN  THE  INDEPENDENT 


VARIABLE  X, SUBJECT  TO  INITIAL  CONDITIONS  YTjT.(J=1.2 NJ. 


CONT.  ON  PO  2 
-lU- 


12 


EflCH  f(J).  THE  QERIVRTIVES  OF  Y(JJ.  MUST  BE  COMPUTED  FOUR  TIMES 


PER  INTEGRATION  STEP  BY  THE  CALLING  PROGRAM.  THE  FUNCTION  MUST 


BE  CALLED' five' TIMES  PER  STEP(PHSSm...PAS6(SJ)  SO  THAT  THE 


n 


INDEPENDENT  VARIABLE  (XJ  AND  THE  SOLUTION  VALUES 


(Y(l) Y(N))  CRN  BE  UPDATED  USING  THE  RUNGE-KUTTR  ALGORITHM 


H IS  THE  PASS  COUNTER.  RUNGE  RETURNS  AS  ITS  VALUE  I -Q  TO 


12 ^ 

SIGNAL  THAT  ALL  DERIVATIVES  (THE  F(JJJ  BE  EVALUATED  OR  0-0  TO 


SIGNAL  THAT  THE  INTEGRATION  PROCESS  FOR  THE  CURRENT  STEP  IS 


FINISHED.  SR'VE'Y(J)  IS  USED  TO  SAVE  THE' 'INITIAL' VALUE  OF  YI'jTI 


AS  arTtt'e'n ."'n  may  be  no  larg'er  than  so.  ) 


CONT.  ON  PO  3 
-115- 


eo.  .^.5. 


CONT.  ON  PO  ♦ 
-IIG- 


eja_3_oF. 


00  50  J=1.3 


1 


1 UPHISC  = -6.2831853071795  Ut^HISC  | 


3 


EC_2_13f_. 


03  TO  200 


COMPUTE  Y COORDINfiTE  OF  SLIPPERS  1 RK'O  4 IH  CCS 


YRSLP  aYCSC«-OSC(  1 .2)kX0555(  1 )<-DSC(  2 .2)kY0SS6(  I )+-DSC(3.2icZ0SS8(  I ) 
YLSLP  =YCSC«-DSC(  1 .2)KXaSS3(43^0SC(2.2J*Y0SSB(4)«-0SC(3.2lBi0SSB(4) 
IF  (YRSLP  .GT.  RRILKT  -fiRD. 


. YLSLP  -GT.  KflILUT)  ISLIP  I 
CRLCULflTE  FORCES  RKO  MOMENTS  RT  ERCH  SLIPPERS 


COMPUTE  POSITION  OF  ITH  SLIPPER  IN  CCS 


XSLIP  iXCSCHlSCC  1 .1  )rXOSSS(  I )•^OSC( 2. 1 IkYOSSBI  I )«-0SC(3.I  )eiOSSf5(  I) 
YSLIP  =YCSCt-nGC(  1 .2)>eXnSSU(  n«-DSC(2,2)KY0SS3(  r)t-DSC(3.2)cI0SS3(  1) 
ISLIP  =ZCSCH>GC(  I .3)bX0SSD(  I )«-DSC(2.3)KYaSS&(  n«-DSC(3.3)B7nSS[J(  I) 


^'cOHPUrl  TUBE  BENDING  fORCES  RNO  HOWENTS 


0B5 


COf.TIKUE 


COflPUIE  PCGITIOM  OP  TUSE  CONfRCT  POINT  IN  CCS 


XCCTCP  =-XCSC  «■  DSC!  I .DkXSTCP  «•  DSC(2.l  )tYSTCP  *■  I’i 
YCCTCP  » YCSC  «•  OSC(  I .E)«XSTCP  DSC( 2.2)kYSTCP  *■ 

ICCTCP  s ICSC  » OSCI 1 .3)»XSTCP  DSG(2«'3)kYSTCP  » 0SC(3«3Ik1STCP 


COMPUTE  TUE£  EXTENSION 


■rL  =.  SQRT(XCCTCF«*2  (VCCTCP  - rUBLTH3)K*2  «•  ZCCTCPi:*2) 
TUSEX I =RL  


RETURN 


r 


I 


CONr.  ON  PO  12 

-10.!)-  EQ-LLJE L"* 


J 


r 


CONT.  ON  PO  3 
-133- 


PQ.2  -OP 1 


0THETR3  = (THETfiS  - THETfiSl )/OELr 


ThETRDl  = TKETRB 


BEKOfiNO  = (UTHESC  - THETflO)a57. 2357795130662 


COMPUTE  MOMENTS  DUE  TO  TDf5S[0.'(RL  SPRINGS 


(YKTUEkCUPCISC 

(PKTUSkC UTHESC 

(RKTUBeCUPHISC 

-YRH 

--PTCH 

-ROLL 


PGISCO)) 

THETfiU)) 

PHtSCO)) 


(YCTUSaOPSI ) 

(PCTU32(0THETfl 

(RCTUbKOPHU 


- OTHETflB)) 


OUTPUT 


V 

K(6.1)  = YflH 
K(6.2)  = PTCH 
11(6, 3J  = ROLL 

' . I T « I. 


12 .V 

COMPUTE  TOTRL  MOMEN  'S 


CONT.  ON  PO 


EQ_3_JJE 


ISRVE  PGSITIOM  OF  SriCG  [ 


7 

SRVECO( IF  INS 
SRVECOI IFINR 
SflVECGdFIKfl 

L.l)  = XCCSM 
L.2)  = YCGSfl 
L.3)  = ICGSri 

8RVE  FORflER  POSITIOKS  RNO  VELGCITIES~i 


< 03  100  1=1.6  > 

7 

I 8flV£Y(  IFIHHL.I)  = Yd)  1 

^ ion  ^ 

1 

1 CONTINUE  1 

12  ^ 

7 

ISRVE  INmfiL  CONDITIONS  | 

\ 

< 00  200 

1=1.121  > 

s , 

1 SflVECQKC  IFINfiL.n  a:  CON(  I ) ] 


COMT.  ON  PO  2 


Pt>  1 OF ?. 


f 

ll  3X.7HIFrNRL=^I3.//)  I 

iRETUf^N  I 


-138- 


PQ  3 FINAL 


CONT.  ON  PO  3 

-140- 


PQ  2 OF { 


loo  TO  520  ^ 


CONT.  ON  PO  ♦ 


1,11- 


CQ.JL_Qf- 


-143- 


Pft  «;  Frw:; 


SUBROUTINE  TRPNSM(QY)  | 


12 i 

L 

1 BEGIN  PDDING  THE  CHPNGE  IN  THE  CENTER  OF  CRRVITY  1 

12  _ Y 

2 

I TO  THE  NECESSRRY  VECTORS  IN  THE  SEPT-HPN  COOROINPTE  1 

12  ^ 

2 

1 SYSTEH.  1 

12  ^ 

TRRNSFORtI  0 

, 

RCG  TO  sues  1 

DKCG-  = D?:CG 
OYCe>  = -DYCG 
DcCO  = -OZCG 
QY116)=  QY(IS) 
QY(17)=  QY(17) 
QY(1Q)=  QY(18) 


PSECl.n^DXCG  - fiSEC2.n»0YCG  - fiSE( 3 . 1 )*DZCG 
PSEC I .ZIsDNCG  - R£ECZ.2)cOYCG  - PSE( 3 .2 JkDZCG 
PSEC  I .S)*DyCG  - RSE(2.3)«:DYCG  - PSE(3.S)yDZCG 


SL 


<rDO  100  I = I .>7  > 


V 

J=3*l-2 

VECTICJ)  rVECTKJ)  - DXCG 
VECTlCJ-fl  )=VECTI(  J+1)  - DYCG 
VECTIC J->2)=VECTU vJ-t-2)  - DZCG 


I CONTINUE  1 


V 


CONT.  on  PG  2 


PG  1 OF  _3- 


-152-  CONT.  ON  PO  2 


Pn  I OF 


-153-  CONT.  ON  PO  3 


EfL2_Q£_--Il 


CONTINUE 


I coNrr'NUE~1 


V 

CONT.  ON  PO  6 

on  C AC  1 


-150- 


CONr.  ON  PO 

Ml  f\- 


10 


EO-S— OE A2 


BLOCltU(  1) 
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APPENDIX  A 


DYNAMIC  RESPONSE  INDEX  (DRT)  EQUATIONS 


DYNAMIC  RESPONSE  INDEX  (DRI)  EQUATIONS 


The  following  spring  clamper  equation  is  utilized  in  caleulating  the  dynamic 
response  index  (DRI),  a measure  of  the  compressibility  of  the  pilot's  spine: 

1.  Y + cY  + kY  = Ay 

Where:  c = damping  coefficient  (/sec) 

k = stiffness  consUinl  (Ibs/in) 

Ay  = vcrlie;il  acceleration  (ins/sec  ^) 

Y = displacement  in  inches 


2.  DRI  = (u)n2/g)  Y 

Where:  =k/m 

m = mass  of  scat /man  combination  (Ibs-sec  ^/in) 

= natural  freq  (rad/sec) 
g = acceleration  due  to  gravity  (in/scc 

The  above  equation  (1)  is  numerically  integrated  and  the  DRI  is  calculated 
using  the  displacement  obtained  at  each  time  step  of  the  ejection  trajectorj’. 


-Al- 


DYNAMIC  CENTKH  OF  GRAVITY  EQUATIONS 


The  displacement  of  the  pilot's  center  of  gravity  from  its  static  position 
been  modelled  using  3 independent  spring-damper  equations,  one  for  each  iixis 
motion. 


has 

of 


1.  Along  the  x-;ixis  of  motion 
X+  c^  X 4 F(x)  = Ax 

I 

Ax  = Fx/M 

Where:  Fx  “ impressed  force  along  the  x-axis  (lbs) 

M = mass  of  seat/man  combination  (Ibs-sec  2/ft) 

Ax=  acceleration  of  scat/man  along  the  x-axis  (ft/sec  2) 

cx  = damping  constant  (/sec) 

F (x)  = elastic  acceleration  (ft/sce  2) 

F (X)  subject  to  the  following  constraints: 

X < 0,  F (X)  = kn  (X) 

0 < X < xslack  , F (x)  = 0 
^ > ^slack»  ^ "^slack^ 

\\3icrc:  x = displacement  along  the  x-axis  (ft) 


Ay  = iicceU'rntion  of  scnt/man  along  (ho  y-axis  (fl/soc 
Cy  = damping  constant  (/see) 

Sy  = spring  modulus  constant  (/see  2) 

3.  Along  the  z-;ixis  of  motion 

Z + C2  Z = F (z)  = A.j. 

Az  = F^/M 

Wliere:  = impressed  force  acting  on  the  seat/man  combination  along  the 

z-:ixis  (lbs) 

M = mass  of  seat/man  (Ibs-sec  ^/ft) 

Az  = acceleration  of  seat/man  ;ilong  the  y-;Lxis  (ft/sec 

Cz  = damping  constant  (/sec) 

F (z)  subject  to  the  following  constraints: 

Z < ZBOT,  F(Z)  = SZN2  * (Z-ZBOT)  + SZNl  * (ZBOT) 

ZBOT  < Z < 0,  F(Z)  = SZNl  * Z 
0<Z<SLACK,  F(Z)  = 0 

ZSLACK  <Z,  F(Z)  = SZP  * (Z-ZSLACK) 

Wlicre  : Z = displacement  along  the  z-axis  (ft) 

ZBOT  = bottoming  distance  in  (he  (-)  Z direction  (ft) 

ZSLACK  = slack  distance  in  () ) Z direction  (ft) 

SZNl,  SZN2,  SZP=  spring  modulus  constants  (/sec  ^ ) 


APPENDIX  C 


MARTIN- BAKER  CATAPULT  EQUATIONS 


L 


MAHTIN-BAKKll  CATAPULT  EQUATIONS 


The  following  linear  spring-damper  and  tube  liending  equations  are  used 
to  simulate  the  Martin-Baker  catapult  system. 

1.  Rail/Slipper  Interaction 

There  are  six  slippers  attached  to  the  scat,  three  on  each  side  of  the 
scat,  starting  at  the  bottom  of  the  scat  (right  hand  side)  - slippers  1, 
2,  and  3,  (on  the  left  hand  side)  - slippers  4,  5,  and  C.  During  the 
catapult  phase  of  the  trajectory,  the  right  and  left  slippers  move  up 
and  off  the  right  and  left  rails  (respectively)  of  the  catapult  (tube  3). 
The  following  spring-damper  equations  arc  used  to  arrive  at  the 
forces  and  moments  acting  on  the  seat/man  during  this  segment  of 
the  catapult  phase: 

^slipj  ^ RAILNT  then  I'sxj»  ^szj*  ^syj*  ^^sxj*  ^^szy»  ^^syj 
otherwise  f'csxi  “ ^s  ^^i  ^^s  ^^i 

^eszj  “ '^s  ^^i  ' *^s  '^^i 
Fcsyi=  l's(Fxi  Fzi 


^SSXj 

F 

CSXj 

^ssYi 

= DSC 

F 

csyj 

^‘sSZj 

F 

^ CSZj 

_ _ 

_ ^ 

Wliere  DSC  is  the  rotation  from  the  catapult  (CCS)  the  seat/man 
coordinate  system  (SMCS) 


— 

■“  — 

— -- 

Fss. 

1 

^sseg. 

ssx. 

1 

^^SSJ 

= 

Y 

sscgj 

X 

^'ssyj 

_Nssi_ 

_ ^SSCgl_ 

SSZj 

-Cl- 


. Torcj  j due  to  the  lorsioniil  sprinfj-diimijcr  egiuilions  ucUng  in  the  jhtch 
plane. 

TOUPTCH  = k'j'OH  A 0 + ‘-'TOH  A 0 

Wlicre:  kg  = slipper  spring  modules  constant  (Ibs/ft) 

Cg  = slipper  damping  constant  (Ibs-sec/ft) 
pg=  slipper  coefficient  of  friction 

Axj  = displacement  of  slipper  (CCS)  from  the  rail  in  the 
x-direction  (ft) 

^Vj  = displacement  of  i^^^  slipper  (CCS)  from  the  rail  in  the 
y-dircction  (upward)  (ft) 

A/.j  = displacement  of  slipper  (CCS)  from  the  rail  in  the 
^--direction  (ft) 

^"csx-*  ^"esv-*  ^'esz-  ~ forces  in  x,  3’,  and  z-direction  (lbs)  in  CCS  due 

to  rail/slipper  interaction 

^ccr.fr  > = distance  from  slippci’;  to  seat/man  center  of 

1 

gravity',  x,  y,  and  z directions  (SMCS)  I’c- 
spectivcly  (ft) 

Lssxi*  ^fssxi*  ^’ssxi  “ moments  produced  by  the  rail/slipper  interaction 

around  the  seat/man  center  of  grav'ity  (SMCS) 
(ft-lbs) 

Utor  = torsional  spring  modules  constant  (Ibs/ft) 
ctor  ~ torsional  damping  constant  (Ibs-scc/ft) 

A 0 = angle  between  r;uls  and  seal  back  (rad) 

AG  = angular  velocity  of  A 0 (rad/scc) 

The  following  tube  extension  equations  :irc  used  to  arrive  at  the  forces  and 
moments  acting  on  the  seat/man  during  the  cataijult  phase  of  the  trajectory. 

Tube  extension  - vector  from  lop  of  tul)e  :l  to  top  of  telescoping  tubes  2 
and  1 (tube  2 contained  inside  tube  .1,  tube  1 contained  inside  tube  2 at  start  of 


-C2- 


1 


calnpuU  phase).  Tube  2 starts  moving  up  and  out  of  tube  3.  Wlien  the  maximum 
length  of  tube  2 is  reached  then  tube  1 starts  moving  up  and  out  of  tube  2. 

TUBEXT  = [Xectep  ^ * (Yeetep  “ TUBLTH:!)  2 4 2]  (CCS) 

If  TUBEXT  > TUBLTlIl  ' TUBL'J'II2  set  all  foi-ces  and  moments  acting  on 
scat/man  due  to  the  catapult  equal  to  /.ero. 

Otherwise,  compute  forces  and  moments  acting  on  seat/man  due  to  the 
tube  bending  as  follows: 

Compute  neflection  of  tube  in  X-Z  plane  (CCS) 

DEFSL=  (Xectep  2 i Zedep  -) 

0 < TUBEXT  < 2 

FTUBE  = ks*DEFSL  + Cs*DEFSL 
OaB  = (I'tUBE  * TITBEXT  2)  / (2Mi*(li  + I2)) 

2 < TUBEXT  < TUBLTUl 

FTUBE  = [3  + DEFSL*E*(Il  + l2)|  / TUBEXT^ 

eAB=  [FTUBE*Tn3EXT  ^]  / [2*E*(I;L  i lo)] 

TniLTH]^<  TUBEXT  <TUBLTll2 

FTUBlf  (>‘^*ItEFLSIUE)  / [(TUBEXT  - TUBUTIlj)  ^ /l^ 

4 (TUBLTlIl  (Ii  t I2))] 

OaB=  |l-rUBE/(2E)]  [ (TUBEXT  - TUBI^TII,)  Vll 

4 TUBLTUl  2 / (i^  4 l2)| 

TUBLTII2  < TUBEXT  < TUBl.THi  ( TUBLTII2 

FtubE  = (2*DEFSIAELM0D)  / I ('I'UBEXT  - TtrBLTIl^) 

+ (TUBLTlIl  Tn3LTIl2  - TUBEXT)  (Ii  • I2) 

I (TUBEXT  - TUBLTUl)  ^/l2  ] 

Oab=  (T-pupp  * I (TUBEXT  - TUBLTl  12)  ^ /ll  • (TUBLTlli  t TUBLTII2- 
TUBEXT)  2/  (ii  . I2)  t ('rUBEXT  - TUBLTlli)  2/I2) 


-C3- 


I'xtubc  = - ^'’lubo  " >^^cctcp  / 
ztube  ~ ~ ^'tube  * ^cctcj)  ^ '^bbSL 
^‘  ytubo  “ ~ ^'^^tubc  * ^tiilic 


sxtube 

xtubc 

sytube 

= DSC 

ztubc 

(CCS SMCS) 

y sztube- 

- ^'ytube- 

^'tube 

^step 

sxtube 

^Hubc 

zz 

^ step 

X 

^sytube 

(SMCS) 

- ^’tube  _ 

_ ^step 

- s ztubc  _ 

Moments  due  to  torsioiiiil  spring-damper  equations  al  the  tube  attachment 
point  acting  in  the  pitch,  roll  and  yaw  directions; 

yaw  = Kytiibc  + + Cytube  * 

ptcli  = Kptiibe  * “ ^ + Cptubc  * ® 

• roll=  Kj.t^,be  * AC'+  C^-tube  * 

Total  forces  due  to  the  catapult: 

G 


xcat  “ 

I 

i=l 

P « T,’' 

ssxj  ‘ sxtube 

G 

ycat  “ 

I 

i=l 

^'ssyj'’'  ^sytube 

G 

zcat  ~ 

I 

^sszj  * ^'s  ztubc 

i=l 


-C4- 


Totnl  moments  duo  to  (lio  oat.ipiiU: 

6 

~ 1 ^■'ssi  ' ^'luho  ' 

i=l 

G 

Meat  = I . M, i PTt:il 

i=l 

G 

Neat  = I N,,.  + . ROLL 

i=l 


Where;  Tl'BKXT  = is  the  veetor  (in  feet)  fi-om  the  top  of  catapult  tube  3 to  the 

top  of  the  telescoping  tubes  2 and  1 (CCS). 

TCBEl/rill,  2 and  3 = lengths,  in  feet,  of  catapult  tubes  ],  2 and  3. 


Nectep’  ^’ec-tep  ^eetep  “ coovdimifcs  (Ct'S)  of  the  extending  tubixs 

2 and  1 in  feet. 

DEFST.  = the  deflection,  in  feet,  of  the  extending  tubes  1 and  2,  from 
the  cataijult  initial  position  in  tlu'  x-/.  plane,  CCS, 

I)EFSL=  rate  of  change  of  DEESI,  in  ft/scc  (CCS). 

FTltBE  = resultant  forc-c  in  llis,  acting  on  the  scat/man  due  to  the 
tube  bending  (CCS). 


UM 


tube 


= tube  bending  friction  coefficient. 


Nxtubc’  and  = conijioncnts  of  the  tube'  bcneling  fore'c 

(lbs)  along  the  x,  y and  z eiircctions  (CCS 

Nsxtubc’  and  = ceunponents  of  the  tube  beneiing  force 

(lbs)  alemg  the  x,  y ;niei  z eiircctions 
(SMCS). 


-Cf)- 


t 


angle  between  tangent  to  the  tube  end  and  vertieal  axis  (Y)  of 
of  fhc  CCS,  radians. 

E = elacticity  modulus  of  lubes  1,  2 and  II,  (Ibs/ft  2). 

TI 

II  (Dio‘^  - 

Dfo,  D20  outside  diameters  of  tubes  1 and  2 

Dfi,  D2i  inside  diameters  of  lubes  1 and  2 
(ft  4). 

^step»  ^'step  '^"step  " coordinates,  in  ft,  or  seat  eonlaet  point  with 

teleseoping  tubes  (SMCS). 

^tube’  ^^tube’  ^’tube  ~ moments  produced  by  the  tube  bending,  around  die 

seat/man  center  of  gra\’ity  (SIMCS)  ft-lbs. 

K , , , K . , , K , , = tube  attachment  torsion;il  spring  modulus  eon- 

ytubc  ptube  rtube 

slants,  ft-lbs/rad,  in  the  yaw,  pitch  and  roll 
di  r e c li  on  s r e s pe  e t i el  y . 

Cytubc*  ^’ptube*  ^rtube  " nttaeliment  torsional  spring-damper  con- 
stants, ft-lbs-see/rad,  iivtlieyaw,  pitch  and 
roll  directions  respectively. 

A*]’*  AO  * A<!>  = angular  displacements  (rad),  of  tlie  seat  back  relative  to 
the  catapult  tube  in  the  yaw,  pitch  and  roll  directions. 

A'l'.  A 0,  A(^  = rate  of  change  of  tlu'  Ai|',  AO,  A'l'  (rad/sec). 


-CO- 


^xc-it*  ^ycnt*  ^’zcat  “ forces  (lbs)  acting  on  the  seat/man  due 

to  the  catapult,  SMCS. 


^-■cat* 


Mc-it,  N’eat  “ moments  (ft-lbs)  produced  by  the  catapult 

around  the  seat/man  center  of  gravity  (SMCS). 


-C7- 


